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Abstract 

We compute the single particle inverse lifetime, evaluated in the supercon- 
ducting state. Within the BCS framework, the calculation can be done non- 
perturbatively, i.e. poles can be found well away from the real axis. We find 
that perturbative calculations are in good agreement with these results, even 
for non-zero impurity scattering. With electron-phonon scattering added to 
the problem, we use the Eliashberg equations, with impurity scattering, to 
calculate the quasiparticle inverse lifetime perturbatively. In all cases we find 
that the inverse lifetime is significantly suppressed in the superconducting 
state, particularly in the presence of impurity scattering. We also compute the 
far-infrared and microwave conductivity, and describe procedures for extract- 
ing scattering rates from conductivity measurements. While these procedures 
lead to scattering rates in qualitative agreement with the inverse lifetime, we 
find that quantitative agreement is lacking, in general. 
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I. INTRODUCTION 

The quasiparticle lifetime is a concept which is useful for clarifying the nature of the 
system of interest, i.e. is it a Fermi vs. marginal Fermi vs. a Luttinger liquid, etc [jl). Various 
techniques are available for measuring lifetimes ; perhaps the most direct is through tunnel 
junction detection ||. In this paper we analyze electron scattering rates, as measured by 
microwave and far-infrared conductivity measurements, and examine their relationship to 
quasiparticle lifetimes. The experimental results for YBaCu307_ x [§]|| and their theoretical 
implications [f)|-§| have been discussed extensively in the literature. We include a treatment 
of inelastic scattering (through the electron-phonon interaction) and use Eliashberg theory to 
further analyze the problem. To some extent this has been done in the normal state already 
PJlO[| . Certain aspects of quasiparticle lifetimes are clarified and the analysis is extended to 
the superconducting state, following our recent work on the optical conductivity [[U||ll[] • 



We will proceed by defining what we mean by a single particle lifetime, bearing in mind 
that sometimes such a definition is not very meaningful, or useful. We review two ways in 
which conductivity data can be used to extract a quasiparticle scattering rate, one based on 
the two-fluid model, and the other based on a straightforward Drude fit to the low frequency 
conductivity. The instances in which these procedures yield a qualitative or quantitative 
facsimile of the quasiparticle inverse lifetime will be clarified. In this way, one can evaluate 
the usefulness of the two-fluid hypothesis, for example, in systems with both elastic and 
inelastic scattering channels. 

This paper has been divided into two sections. The first reviews quasiparticle lifetimes, 
which are calculated from the single particle Green function. Previous work is extended to 
include both electron-impurity scattering (in the Born approximation) and electron-phonon 
scattering. Perturbative methods within the Eliashberg framework are used. The second 
section elucidates two methods recently used to extract scattering rates from optical and 
microwave conductivity measurements, and how these relate to the single particle inverse 
lifetime. We should make it clear at the start that the calculation we use for the opti- 
cal conductivity omits vertex corrections. The main result of this is that the "1 — cos 9" 



factor that occurs in a Boltzmann formulation of transport properties is absent fll2| |, so in 



fact, when theoretical comparisons between the quasiparticle lifetime and the conductivity- 
derived scattering rate are made, the agreement will in general be better than it would be, 
given an exact calculation. This is made obvious in the next section, which treats elastic 
impurity scattering. Conversely, good agreement between these two quantities in the theory 
does not imply good agreement between the two measured quantities. 
II. QUASIPARTICLE LIFETIMES 

II. A Normal State 

The fact that the quasiparticle lifetime is related to the conductivity is most apparent 
in the Drude model, which incorporates electron-impurity scattering into the problem. In 
the Born approximation, the single-electron self-energy is given by 

E(a; + i5 ) = _J_ (1) 

where 1/r is the impurity scattering rate and r is the single particle lifetime. We will utilize 
an expression for the optical conductivity which omits vertex corrections, and so it is given 



simply as the convolution of two single particle Green functions fl"3|-|T6f . The result reduces 
to the Drude conductivity [|T7j], 

a{u) = — 2 

m 1/r — iv 

whose real part is given by the usual Lorentzian form with half- width 1/r. Note that vertex 
corrections would actually alter the 1/r that appears in eq. (|2|); however, if only s-wave 
phase shifts are included, the two inverse lifetimes would be identical. Here n, e, and m 
are the free electron density, charge magnitude, and mass. Thus the single particle lifetime 
appears explicitly (and can easily be measured) in the optical conductivity. 

The next step is to include phonons, and thus an additional scattering mechanism, via 
the electron-phonon interaction fll8|,|9|| . Then the electron self-energy is given by fl9|j 
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(3) 



where a 2 F(v) is the electron-phonon spectral function, N{y) is the bose function, and i/j(x) 
is the digamma function. The quasiparticle has energy and lifetime defined by the pole of 
the single particle retarded Green function, i.e. the zero of 

G~ l {k,uj + i8) = u; - e fc - + i<5), (4) 

where now the functions are analytically continued to the lower half-plane. Defining the 
energy and inverse lifetime by 

u = E k - iT k (5) 

we obtain the two equations, 

y°° 2 { I v-Ek+J^k, ,A .v + E h -iT k \ 
E k = e k + J^ dva F{u)Rei^{- + i — ^ )> (6) 

r fe = -I + r du^F(J2n[N(u) + V ImUi 1 - + _ ^ _ i ^E k -iT k , 

k 2r io y ! \ [ y! T l^ v 2 2%T ; ^ y 2 2%T ^ 

(7) 

The general problem requires a numerical solution of these equations and is beyond the 
scope of this paper. Here, we assume that E k and T k are small so that one can obtain closed 
expressions. They are [ |19| . ^0[ 



Ek ~ i + aV) (8) 

k ~ 1 + A*(T) [) 
where f(v) is the Fermi function and 

A*(T) = -— jf dvc?F{v)Im^{~ + i— ). (10) 

As T -> 0, A*(T) -> A = 2 / °° ^a 2 F(z/), while at high temperature, A*(T) oc 1/T 2 . Typical 
results for various electron-phonon spectral functions are shown in Fig. 1. 



II. B Superconducting State 

Quasiparticle lifetimes in electron-phonon superconductors were discussed long ago by 
Kaplan et al. 0. Before proceeding to their relationship with the optical conductivity we 
give a brief review of the formalism 0,2~I]. As in the normal state the quasiparticle scattering 



rate is defined by (twice) the imaginary part of the pole in the single particle Green function. 
In the superconducting state, the diagonal component of the single particle Green function 
2| 



is 



where Z(u) and 4>{oj) are the renormalization and pairing functions given by solutions to 



the Eliashberg equations [23,22.24] which are repeated here for convenience: 
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and 



uo{uj) = uo + nrl ~ luJ n 
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o I ' ^Cj 2 (uj -u)-<p 2 (uj- v) 

^Juj 2 (uj + v)- (j) 2 (uj + v) J 
where u(u) = uZ(u). Here, N(y) and f(u) are the Bose and Fermi distribution functions, 
respectively. The electron-phonon spectral function is given by a 2 F(u), its Hilbert trans- 
form is X(z). The Coulomb repulsion parameter is fi*(uj c ) with cutoff u c . A negative value 
for this parameter can be used to model some BCS attraction of unspecified origin. The 
renormalization and pairing functions are first obtained on the imaginary axis at the Mat- 
subara frequencies, i.e. uo = iuo n = mT(2n — 1), with <p m = <p(iuj m ) by setting the complex 



variable uo in these equations to the Matsubara frequencies p5[^6j. Then the equations are 



iterated as written, with u set to a frequency on the real axis. Note that the square roots 
with complex arguments are defined to have a positive imaginary part. 

Actually, these functions are obtained at frequencies just above the real axis, and then, in 
principle, the pole is given by the zero of the denominator continued to the lower half-plane. 
However, since the solutions are readily known only along certain lines in the complex plane 
(e.g. the imaginary axis or just above the real axis) we follow previous authors and look for 
the pole perturbatively. That is, we write u = E — iT, and linearize the imaginary part so 
that MM 



T(E) = ~UE) EZf(E) (M) 

Note that E is determined by equating the real parts (again linearizing in imaginary com- 
ponents), which yields 



E 



(15) 



At the Fermi surface, e^ = so E — <fix(E) / Z\(E) = Ai(E) becomes the definition for the 
lowest energy excitation, i.e. the gap in the excitation spectrum. It has become common 
practice to consider E as an independent variable, and then to study T(E) as a function of E. 
In Fig. 2 we show T(E) vs. E for various temperatures, using a Debye model spectrum. It is 
clear that the scattering rate near the gap edge (shown by the arrow) decreases very quickly 
as the temperature decreases. Note that T(E) actually becomes negative at intermediate 
temperatures near E « 10 meV. This is real (i.e. not a numerical artifact) and is simply a 
property of the function T(E) given by eq. ([[4]) through the linearization procedure. The 



true pole in eq. ( |iT|) will always have a negative imaginary part (i.e. T is positive). 

What happens in the BCS limit ? Then Z 2 (E) — > and <p2{E) — > 0, i.e. the functions 
involved in the solution are pure real, with Z\{E) — ► 1 and <fii(E) — ► A(T), where A(T) is 
obtained self-consistently from the BCS equation. Clearly then the quasiparticle energy is 



E = J el + A 2 and the scattering rate, T = 0. Thus, within the BCS approximation the 



quasiparticle states are infinitely long-lived, as there is no means by which a quasiparticle 
can decay. In the electron-phonon theory, the phonons provide the means. In this way the 
BCS approximation is pathological. 

Before including the electron-phonon interaction as a scattering mechanism, one can first 
introduce impurities to provide elastic electron scattering. Then 

4>(UJ) = 4> c l{u) + 



2t 



^UJ 2 Z 2 (u) - 2 (u) 



Z{u) = Z cl {u) + 1- Z{UJ) (16) 

2t ^lo 2 Z 2 (u) - <P 2 {uj) 

where the subscript £ cl' refers to the clean limit and 1/r is the (normal) impurity scattering 
rate. The gap function is defined 

A(w) = 4>(u>)/Z(u) (17) 

and is independent of impurities. In the BCS limit, 

i Asgnu; 



(f>(u) = A + 



2r Vcj 2 - A 
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Z(u) = l + ^-^E^. (18) 
V ; 2r Vcu 2 - A 2 V ; 



As before, the quasiparticle energy is given by E = \J t\ + A 2 , but the scattering rate is now 

(19) 
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2r \E\ 2r^ + Aa 

Eq. (|T9|) shows that at the Fermi surface impurities are ineffectual for electron scattering, 
i.e. this is another manifestation of Anderson's theorem This happens abruptly at T c ; 



as soon as a gap develops the scattering rate becomes zero. Away from the Fermi surface the 
scattering rate is reduced from its value in the normal state. In particular, well away from 
the Fermi surface the scattering rate in the superconducting state approaches the normal 
state rate. This behaviour is illustrated in Fig. 3. 

The calculations in Fig. 3 rely on a perturbative search for the pole in the lower half- 
plane, as given by eq. (|14|) . A more rigorous calculation reveals poles "within the gap", but 



these do not appear in the spectral function because of coherence factors. This can be seen 
by rewriting eq. ([LJ]) as 

G u (k,u) = Ul JJ X * 1 



2 V V^ 2 - A 2 J Vuj 2 - A 2 - e k + ^sgna; 
Ml- , " ) , — (20) 

The coherence factors, ± J , are such that the spectral function, A(k, u) = 

—-ImG(k,uj), is zero for |a>| < A, independent of the impurity scattering rate, as can be 
readily verified explicitly from eq. (pUp . On the other hand the pole of the Green function 
is given by solving for u = E — iT in the two coupled equations 

^ 2 -r 2 = A 2 + 4-(^) 2 (2i) 

It 



(22) 





1 


\E\ 


2r 


Eh 





'e\ + A 2 , and V — 0. For finite impurity 
scattering, however, the solution depends on at the Fermi surface (e& = 0): 



E = 4 /A 2 - (— ) 2 , T = for M- < 1 (23) 

V 2r ; ' 2rA v ; 



E = 0, r = ^) 2 -A 2 for ^a>1 (24) 

In either case the solution is irrelevant, because the real part lies within the gap (i.e. the 
quasiparticle residue is zero). The relevant energies are \E\ > A, and then the scattering 
rate given by Eq. (|22"D agrees with the linearized solution, Eq. (|14"D, although the quantity V 
no longer truly represents a quasiparticle inverse lifetime. Away from the Fermi surface the 
quasiparticle energy requires the solution of a quadratic equation and E may lie below or 
above the gap, depending on e k and 1/r. Once more solutions with \E\ < A are irrelevant 
because the residue is zero. In Fig. 4a (b) we show the real (imaginary) parts of the pole for 
some typical parameters. As can be seen, some atypical behaviour can occur as a function 
of temperature, for example when e k = 0. As the temperature is lowered the pole moves 
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from a point on the negative imaginary axis to the origin (near T/T c = 0.8) and then moves 



along the real axis towards E — y A 2 — (^:) 2 . Nonetheless, comparison of Fig. 4(b) with 
Fig. 3 shows that the perturbative is in qualitative agreement with the non-perturbative 
solution. 

Coherence factors have always been thought to play an important role for two-particle 
response functions (such as the NMR relaxation rate or the microwave conductivity). In 



particular they lead to the Hebel-Slichter p8| singularity in the NMR relaxation rate. How- 



ever they are considered not to play a big role in single particle functions. Eq. ( |20"D shows, 
however, that the coherence factors play a very important role in the single particle spectral 
function, in that they maintain a gap equal to A, even when impurities cause the single 
particle pole to have a real part whose value lies in the gap. 

The frequency dependence of the spectral function is shown in Fig. 5. For energies close 
to the Fermi surface the spectral function is dominated by the square root singularity at the 
gap edge, which comes from the coherence factors rather than the single particle pole. For 
larger energies the peak present is due to the pole and differs from the normal state spectral 
function only at low frequencies. 

We return now to the strong coupling case. In the superconducting state, we do not have 
access to the analytic continuation of the Green function to the lower half-plane p9| . We 
thus follow Karakozov et al. ||30||, and use an expansion near u = 0. We also follow Kaplan 



et al. H and henceforth use Eq. ([14]) for the scattering rate, noting that while the energy 
E at which T(E) is evaluated ought, in principle, to be computed self-consistently as was 
done in BCS, i.e. eqs. ( [21~| , |22"|) , in practice we will choose the relevant energy, and so treat 
energy as an independent variable. 

Karakozov et al. |30| pointed out that at any finite temperature the solutions to the 
Eliashberg equations have the following low frequency behaviour: 

Z{u) « + ^ (25) 
to 

A(u) « <V 2 - i5 2 uj (26) 
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where Z"i, 72,^1 and b% are real (positive) constants. This implies that a quasiparticle pole 
exists with 

E « ^= (27) 



r « ^(2vr r dva 2 F{v)g{v)[N{v) + /(v)] + ^) (28) 
Zi v JO It ' 

where we have used Eq. (|14]) and assumed E — > (i.e. we are near the Fermi surface). In 
Eq. (|28"D g(z/) is the single electron density of states in the superconducting state, 

jM = i K7ra)' (29) 



which is non-zero for v = at any finite temperature [p0| -|32]|. In Fig. 6 we show the 
quasiparticle scattering rate, 2T(T,E = 0) vs. T/T c for various impurity scattering rates, 
1/t, in the superconducting state. The normal state result is also shown for reference. Note 
that in the clean limit there is an enhancement just below T c in the superconducting state, 
but at low temperatures there is an exponential suppression, compared to the power law 
behaviour observed in the normal state. When impurity scattering is present there is an 
immediate suppression below T c (a knee is still present, however, in the superconducting 
state). It is clear that at low temperatures the superconducting state is impervious to 
impurity scattering, as one would expect. 

Is 2r(T, E — 0) the relevant inverse lifetime by which properties of the superconducting 
state can be understood ? The answer is no, at least not at low temperatures, where the 
spectral function essentially develops a gap for low energies. In Fig. 7 we plot the spectral 



function at the Fermi surface, A(kp,ui) for various temperatures | 3~1|| . While no true gap 
exists at any finite temperature, it is clear that for low temperatures the spectral weight at 
low frequency is exponentially suppressed. Thus we have a situation similar to that in BCS 
theory, where the imaginary part of the pole is irrelevant; rather it is the scattering rate 
(as defined by Eq. flUD) evaluated above the "gap-edge" where a large spectral weight is 
present, that is most relevant. When impurities are added the spectral peak is broadened, 
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even at low temperatures, as shown in Fig. 8 for an intermediate temperature, T/T c = 0.5. 
However the lineshape is very asymmetric as a gap remains at low frequencies (particularly 
prominent at low temperatures). Fig. 7 demonstrates that the relevant energy is a function 
of temperature. In Fig. 9 we plot 2T(T,E = A(T)) vs. T/T c , where A(T) is determined 
from the relation (at the Fermi surface) 

E = ReA(E, T). (30) 

By Eq. (|3UJ) we understand that the E — solution is excluded (similarly the very low 
energy solution is also excluded). We are interested in the conventional solution which gives 
rise to the peak in the spectral function illustrated in Fig. 7 or 8. If no non-zero solution 
is present (as is the case near T c ) then we utilize the E = solution. Note that in this 
case care must be taken when obtaining the E — > limit of Eq. (|14]) . Also there is present 
in this definition a discontinuity at some temperature near T c , which is where a nonzero 
solution first appears. In this way we hope to show the scattering rate at an energy where 
the spectral weight is large, and therefore of most relevance to observables. At any rate, it 
is clear by comparing Fig. 9 to Fig. 6 (see also Fig. 2) that there is very little difference in 
the scattering rate in the gap region of energy. However, as the energy increases beyond the 
gap the scattering rate increases, resulting in short "lifetimes" , even at zero temperature || . 

III. THE OPTICAL CONDUCTIVITY 



Using the Kubo formalism |12| , the optical conductivity can be related to a current- 
current correlation function. The final result for the frequency dependence of the conduc- 



tivity in the long wavelength limit is |jl3|-|16 



<r(u) = -(u(v + iS)+—), (31) 

where the paramagnetic response function, n(z/ + iS), is given by 

H(u + iS) = 1 — 1 + / du ta.nh( — )[hi(uj,u + v) — h 2 (uj,u + v) ) 

m [ Jo 2 v ' 

+ J° duo tanh( ^ 2 + U ^ )(h\{u, uo + v) + h 2 (u , uj + i/)) j (32) 
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with 



1 - NiuJNfa) - PMPfa) 

ftl(Wl,W 2 ) = — — — TT 

2(e(cJi) + e{uj 2 )) 
1 + N'MNfa) + P*(cui)P(cu 2 ) 



^2(^1,^2) 
iV(u;) 



2(e(w a ) - e-M) 



e(w) = ^Cj 2 (uj + i5)-<p 2 (uj + i5). (33) 

In eq. (|3~2"D D is a large cutoff of order the electronic bandwidth. 

To summarize the previous set of equations: for a given model for the spectral density, 
a 2 F{v) and a choice of impurity scattering rate 1/r, we can compute the conductivity <j(v) at 
any frequency and temperature, including effects due to both elastic and inelastic scattering 
mechanisms, the latter being determined by the choice of electron-phonon spectral density. 
As explained in the introduction, however, vertex corrections are omitted. 

In the analysis of experimental data it is possible to use several methods to extract a 
single temperature dependent scattering rate. One such method utilized the microwave con- 
ductivity in YBaCu307_ x @,|] and in Nb [0], and adopted a two-fluid model description 
of the superconducting state. It was assumed that the absorptive component of the con- 
ductivity (real part of a at finite frequency, denoted by <Ti(u)) was due only to the normal 
component of the fluid. Thus an expression of the form 



■2 



_,, ne m 

°\ \y, T) 



m m*(T) 



1 - 



A 2 (0) 



T 



(T) 



1 + (ur(T)y 



(34) 



A 2 (T)_ 

is assumed to hold approximately. In eq. ([34]) r(T) has units of a scattering time and A(T) is 
the penetration depth at temperature T. Here we will calculate <J\{y, T) for a model a 2 F(v) 
and 1/r using the full expression (|3l]). At the same time we can calculate the penetration 
depth, either from a zero frequency limit of eq. ([31]), or directly from the imaginary axis 
EH: 
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The idea is to examine the zero frequency limit of eq. (|34|) , and thus define a scattering rate 
relative to the rate at T c J/J: 

r(T c ) _ ( A 2 (0) \ a N (T c ) 

r(T) " \ l \\T)j a^T) ' W 
where it has been assumed that the mass enhancement factor in eq. does not change 
with temperature. 

In Fig. 10 we show results for t(T c )/t(T) defined by eq. (^) in the superconducting 
state (solid curve) with which we compare the inverse lifetime for both the normal (short- 
dashed curve, eq. @ ), and the superconducting (long-dashed curve, eq. fll~4] ) ) states. 



The results are based on a Debye model spectrum for a 2 F(z/) used in the previous section, 

with mass enhancement parameter A = 1 and T c = 100 K. (A negative fi* is required.) 

Note that the normal state scattering rate (short-dashed curve) is only sensitive to the low 

frequency part of a 2 F(u) at low temperatures: a v 2 dependence in a 2 F(u) implies a T 3 

dependence in l/r(T) (Again, vertex corrections would alter this to the familiar T 5 law 

||17|| .) The agreement between the inverse lifetime (long-dashed curve) and the scattering 

rate defined by the two-fluid model (solid curve) in the superconducting state is remarkable. 

This indicates that the two- fluid description makes sense and the quantities shown in 

Fig. 10 apply to the normal component of the superfluid. Fig. 10 illustrates the comparison 

in the clean limit, where the two-fluid description is expected to be most accurate ||. Before 

investigating impurity dependence, we turn to a second possible procedure for extracting 

a scattering rate from conductivity data |3~^ , |35|| , which is simply a generalization of that 

used by Shulga et al. [[J to the superconducting state. One simply fits the low frequency 

absorptive part of the conductivity to a Drude form: 

/ m\ n e 2 1 r*(T) 

ai(i/,T) = — (37) 

m m* jm 1 + (z/r*J z 

As described in Refs. [p|JT0[1 , it is possible to fit a Drude form to the low frequency part of the 
optical conductivity in the normal state. Such a fit is also possible in the superconducting 
state [^]. Theoretically, the fit is problematic in a BCS approach because there is a Hebel- 
Slichter logarithmic singularity at low frequency at all temperatures in the superconducting 
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state. However, with the Eliashberg approach, the Hebel-Slichter singularity is smeared, and 
one can fit a Drude form over a limited range of frequency. Such a fit for a Debye spectrum 
is also included in Fig. 10 (dotted curve). While the fit is in qualitative agreement with the 
inverse lifetime, this method of characterizing the inverse lifetime in the superconducting 
state is clearly not as accurate as the two-fluid model. The fits themselves are shown in Fig. 
11. It is clear that the fits fail at sufficiently high frequency, as one would expect, but that 
they characterize well the low frequency response in the superconducting state. 

Fig. 10 clearly shows that there is considerable freedom and hence ambiguity in extract- 
ing a temperature dependent scattering rate from conductivity data. Nevertheless, in the 
clean limit it is evident that the two-fluid model is a useful device for extracting the low 
energy quasiparticle inverse lifetime in the superconducting state. 

With the addition of impurities the situation changes considerably. This is illustrated in 
Fig. 12, where the same calculations as in Fig. 10 are shown, but with an additional impurity 
scattering, 1/r = 25 meV included. The use of formula fl34"|), inspired by the two- fluid model, 
gives a scattering rate (solid curve) that falls much less rapidly around T = T c then does the 
inverse lifetime in the superconducting state (long-dashed curve). The latter curve drops 
almost vertically as the temperature drops below T c , as has already been discussed. The 
result based on the Drude fit (dotted curve) shows a peak which is reduced in size from that 
in the clean limit (Fig. 10). Indeed, for increased impurity scattering, the peak just below 
T c disappears. In both rapid suppression is expected just below T c : in the case of 

the inverse lifetime, this suppression is a consequence of Anderson's theorem | 2"7fl , as already 



discussed. In the case of the Drude fit, it is easy to see that this is the case in the dirty limit. 
In the dirty limit the conductivity is almost flat as a function of frequency on the scale of 
1/r, at T c . Just below T c , however, a gap in the spectrum begins to develop, so that weight 
is shifted from roughly the gap region to the delta function at the origin. Any low energy 
fit will then use a Lorentzian width which senses this depression in the conductivity, which 
is on an energy scale of the gap. This represents a significant suppression from the normal 
state scattering rate (infinite in the dirty limit). Here we are in an intermediate regime, 
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with 1/r = 25 meV (note: A(T = 0) = 20.2 meV). The corresponding fits for Fig. 12 are 
shown in Fig. 13. 

Finally, we examine the low frequency conductivity vs. reduced temperature, a quantity 
measured in the high T c oxides by microwave techniques. In Fig. 14 we show the real part 
of the conductivity, cri(z/) vs. reduced temperature, for several low frequencies. Note the 
relative insensitivity to frequency, a feature of strong coupling pointed out in Ref. |36]. Also 



note the lack of a coherence peak just below T c . Nonetheless, a broad peak exists at lower 
temperatures, somewhat reminiscent of that observed in YBaCu^Oj^x [[37],|]|§. This peak 
exists because of a competition between an increasing scattering time (making a increase) 
and a decreasing normal component (making a decrease, particularly as T — > 0). We should 
note that this peak is most prominent in the clean limit. As Fig. 13 indicates (see values at 
the intercept), the peak is absent for a sufficiently large impurity scattering rate. 

It is of interest to examine what dependence these results have on the electron-phonon 
spectral function. As an extreme we utilize a spectrum which is sharply peaked at some high 
frequency, and, in contrast to the Debye model employed above, coupling to low frequency 
modes is absent; such a spectrum models a strong coupling to an optic mode. We choose a 
triangular shape for convenience, starting at u = 34.8 meV with a cutoff at uje = 35.5 meV. 
The coupling constant is chosen so that the mass enhancement value is A = 1, in agreement 
with that chosen for the Debye spectrum. As was the case there, a negative /i* is used to 
give T c = 100 K, and the zero temperature gap was found to be close to the Debye value. 

In Fig. 15 we plot the real part of the low frequency conductivity, cr 1 (u) vs. reduced 
temperature, now using the triangular spectrum for a 2 F(i>). In contrast to the result for 
the Debye model, the low frequency conductivity is strongly frequency dependent at low 
temperatures. In fact, for v = 0, the conductivity appears to diverge. We believe that 
at sufficiently low temperature, this curve will actually achieve a maximum and approach 
zero at zero temperature, but we have been unable to obtain this result numerically. Once 
again there is a competition between an increasing scattering time and a decreasing normal 
density component as the temperature is lowered from T c . Here, however, the increasing 
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scattering time appears to be overwhelming the decrease in normal fluid density. The key 
difference with the Debye spectrum is that here the spectrum has a big gap, so that the 
lifetime is increasing exponentially with decreasing temperature, already in the normal state. 
Recall that in the Debye case the increase followed a power law behaviour with decreasing 
temperature. Since the decrease in normal fluid density is always exponential, this term 
dominates in the Debye case, whereas in the case of the gapped spectrum the competition is 
subtle, and will depend strongly on the details of the spectrum (an electron-phonon spectrum 
with a much smaller gap will yield a zero frequency conductivity which approaches zero at 
zero temperature, for example). 

The physics of this conductivity peak is different from what has been already proposed 
for YBaC^Oy-a;. It has been suggested that the excitation spectrum becomes gapped due 
to the superconductivity, i.e. a feedback mechanism exists which creates a low frequency 
gap in the excitation spectrum as the superconducting order parameter opens up below T c . 
Such a scenario has been explored within a marginal Fermi liquid scheme [p7|-|39| , and is also 
consistent with thermal conductivity experiments [iOj. Here the conductivity peak arises 



because the spectrum is already gapped, and the scattering rate is sufficiently high at T c 
because T c itself is fairly high. We should warn the reader that this mechanism requires a 
"fine tuning" of the spectrum, i.e. a large gap is required. 

Note that for any non-zero frequency the conductivity has a visible maximum, and 
quickly approaches zero at sufficiently low temperature. Nonetheless this turn around occurs 
for yet another reason: as one lowers the temperature the Drude-like peak at low frequencies 
gets narrower while at the same time the magnitude of the zero temperature intercept 
increases. For any given finite frequency, then, a temperature is eventually reached below 
which this frequency is now on the tail of the Drude-like peak. This means that while the 
zero frequency conductivity increases, that at any finite frequency will eventually decrease 
as the width becomes smaller than the frequency. So in this case the increase in scattering 
lifetime still dominates the decrease in normal fluid density, but, because we are fixed at a 
finite frequency, the conductivity decreases. Note that frequencies of order 0.01 meV (~ 2.4 
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GHz) are within the range of microwave frequencies that are used in experiments. 

To show this more explicitly we illustrate in Fig. 16 the various scattering rates obtained 
with the gapped spectrum. These are to be compared with those shown for the Debye model 
in Fig. 10. Clearly the quantitative agreement between the scattering rate inspired by the 
two-fluid model and the inverse lifetime that was seen in Fig. 10 with the Debye spectrum 
was fortuitous. While a qualitative correspondence between these two entities continues to 
exist with the gapped spectrum, they are no longer in quantitative agreement. We have 
verified that this is generically true, by investigating other spectra, not shown here. 
IV. SUMMARY 

We have investigated the quasiparticle lifetime in an Eliashberg s-wave superconductor, 
generalizing earlier work || to include impurity scattering as well. We find that the quasipar- 
ticle lifetime becomes infinite at low temperatures, independent of the impurity scattering 



rate, which we understand as simply a manifestation of Anderson's theorem [27]. Thus, on 
general grounds, within a BCS framework, the scattering rate should collapse to zero in the 
superconducting state. 

We have also investigated two methods for extracting the scattering rate from the low fre- 
quency conductivity. One relies on a two-fluid model picture, and the other simply utilizes a 
low frequency Drude fit. Neither should necessarily correspond very closely to the quasipar- 
ticle inverse lifetime, and we find that in general they do not, quantitatively. Qualitatively, 
however, the scattering rate defined by either procedure gives the correct temperature de- 
pendence for the inverse lifetime. In the presence of impurities, the two-fluid prescription 
appears to be less accurate, presumably because such a prescription takes into account only 
the lower normal fluid density as the temperature is lowered, and not the fact that impurity 
scattering is less effective in the superconducting state. Thus we caution that interpreting 
a conductivity-derived scattering rate as a quasiparticle inverse lifetime can lead to inaccu- 
racies. 
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ure Captions 



Fig. 1. (a) Mass enhancement parameter, as defined by Eq. (|T0| ) in the text, vs. 
reduced temperature T/T c , for various electron-phonon spectral functions. In all cases 
T c = 100 K. The 4 spectra used are a Debye spectrum (solid line), linear spectrum 
(dotted line), a spectrum proportional to \fv (dashed line), and a triangular spectrum 
(dot-dashed line). In all cases the strength is such that A = 1. In the first three cases 
a cutoff frequency equal to 30 meV was used. In the last case, the spectrum starts at 
34.8 meV and is cut off at 35.5 meV. 

(b) Inverse lifetime, 2Y(T) (in meV) vs. reduced temperature for the same spectra as 
in (a). Note the different low temperature behaviour, depending on the low frequency 
characteristics of the electron-phonon spectrum used. 

Fig. 2. Scattering rate, 2T(T, E) (in meV) vs. E (in meV) for various temperatures in 
the superconducting state. The zero temperature gap at the Fermi surface is indicated 
by the arrow. Note that for T/T c = 0.5 the function plotted actually becomes negative 
(near 10 meV). The physical pole occurs at an energy, E, where 2T(T,E) is always 
positive, however. 

Fig. 3. Scattering rate normalized to the zero temperature gap, 2T(T,E)/A vs. 
reduced temperature T/T c , for various quasiparticle energies, e. These are computed 
using the perturbative expansion, Eq. (|i"9|). The impurity scattering rate is l/(r)A = 
1. Note that on the Fermi surface (solid curve), the scattering rate is zero immediately 
below T c . 

Fig. 4. (a) Real and (b) imaginary parts of the quasiparticle pole vs. reduced temper- 
ature, T/T c , with 1/(t)A = 1 and various quasiparticle energies. These are computed 
non-perturbatively from Eqs. (|21|j22] ). Note that at the Fermi surface the quasiparticle 



energy has an abrupt onset at a temperature somewhat below T c . The results in (b) are 
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similar to those in Fig. 3, except that at the Fermi surface, the decrease in scattering 
rate below T c is not as abrupt. 

Fig. 5. The single particle spectral function, A(k,u), vs. normalized frequency, uj/ A, 
with 1/(t)A = 1, for various quasiparticle energies, e/A. The Fermi surface result 
(solid curve) has particle-hole symmetry. The result for e/A = 0.5 would have a peak 
within the gap (between -1 and 1) except that the coherence factors in Eq. ( ^0|) give 
zero residue for the gap region. As the quasiparticle energy increases, the spectral 
function begins to resemble the normal state spectral function. Small square-root 
singularities still exist, nonetheless at the particle and hole gap edges. 

Fig. 6. The quasiparticle scattering rate, 2T(T,E = 0) vs. T/T c for various impurity 
scattering rates, 1/r, in both the superconducting and normal (long-dashed curves) 
states. In the clean limit (solid curve) there is an enhancement immediately below 
T c . When impurity scattering is present, this scattering is immediately suppressed 
in the superconducting state, as shown by the dotted and dashed curves. Below 
a temperature of about 0.8T C the scattering rate is independent of the amount of 
impurity scattering present in the normal state. A Debye electron-phonon spectrum 
was used with A = 1 and cutoff frequency, ojd = 30 meV. 

Fig. 7. The single particle spectral function, A(kp,uj) at the Fermi surface, vs. fre- 
quency, in the clean limit, for various reduced temperatures. The spectral function is 
considerably broadened near T c , due to temperature alone. At the two lowest temper- 
atures shown, while there is no true gap in the excitation spectrum, this plot makes it 
clear that, practically speaking, an effective gap in the excitation spectrum is present. 
A Debye electron-phonon spectrum was used with A = 1 and cutoff frequency, up = 30 
meV. 

Fig. 8. The single particle spectral function, u) at the Fermi surface, vs. fre- 

quency, for various impurity scattering rates, as indicated. Note the broadening which 
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occurs with increasing impurity scattering. However, spectral weight remains absent 
in the "gap region". Results are for a temperature T/T c = 0.5, and with a Debye 
electron-phonon spectrum with A = 1 and cutoff frequency, ujd = 30 meV. 

Fig. 9. The quasiparticle scattering rate, 2T(T, E = A(E)), evaluated at the quasipar- 
ticle energy, given on the Fermi surface by E = A(E), vs. T/T c for various impurity 
scattering rates, 1/r, in both the superconducting and normal (long-dashed curves) 
states. These results are in quantitative agreement with those in Fig. 6, since the en- 
ergy scale, A, is still small compared to other (phonon) energy scales in the problem. 
In particular, below about 0.8T C , the scattering rate is independent of the amount of 
impurity scattering present in the normal state. A Debye electron-phonon spectrum 
was used with A = 1 and cutoff frequency, ujd = 30 meV. 

Fig. 10. Various normalized scattering rates vs. reduced temperature. The normal 
and superconducting scattering rates come from the quasiparticle inverse lifetime, 
given by eq. fll4]) , at zero energy. The two- fluid result comes from eq. (|36| ) while 
the Drude fit is obtained by fitting eq. (|37| ) to the low frequency conductivity in the 
superconducting state. Note the agreement of the scattering rate as extracted from 
the two-fluid analysis with the inverse lifetime in the superconducting state. A Debye 
electron-phonon spectrum was used with A = 1 and cutoff frequency, ojd = 30 meV. 

Fig. 11. The low frequency conductivity in the superconducting state (solid curves) 
along with their fits based on eq. fl37|) (dashed curves). These fits were used in Fig. 
10 (dotted curves). 

Fig. 12. Same as Fig. 10, except now with an impurity scattering rate of 25 meV. 
Note that the two-fluid analysis agrees poorly with the quasiparticle inverse lifetime. 

Fig. 13. The low frequency conductivity fits used in Fig. 12. 

Fig. 14. The very low frequency conductivity as a function of reduced temperature, 
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in the clean limit. A Debye electron-phonon spectrum was used with A = 1 and cutoff 
frequency, wo = 30 meV. Note that the results are relatively insensitive to frequency 
(the v = 0.01 meV result, given by the dotted curve, is essentially hidden by the zero 
frequency result). A microwave experiment yields essentially zero frequency results. 

• Fig. 15. Same as for Fig. 14, but now with the triangular spectrum as described in 
the text. Note that the zero frequency conductivity appears to diverge as T — > 0. A 
microwave experiment will yield a large low temperature peak as a function of reduced 
temperature, whose magnitude will depend strongly on the frequency. These results 
are for the clean limit. The peak is also reduced as impurity scattering is added (not 
shown) . 

• Fig. 16. Comparison of scattering rates (as in Fig. 10) calculated for the triangular 
spectrum, in the clean limit. Note that the normal state result approaches T = 
exponentially due to the gap in the o?F{y) spectrum (in Fig. 10 the corresponding 
curve approached zero with a power law behaviour). While the results in the super- 
conducting state are qualitatively similar, they no longer agree quantitatively with one 
another, as in Fig. 10. 
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